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Torsion of Anisotropic Elastic Cylinders by 
Forces Applied on the Lateral Surface^ 

Harold Luxenberg^ 

The classical Saint- Venant theory of torsion presents a method for the determination 
of the elastic behavior of an orthotropic elastic beam of uniform arbitrary cross section 
twisted by forces applied on the end sections, with the lateral surface free of stress. In the 
present paper the theory is extended to include a more generally anisotropic beam (possessing 
only a single plane of elastic symmetry, rather than three mutually orthogonal planes of 
elastic symmetry, as in the case of orthotropic materials) built in at one end and twisted by 
forces uniformly distributed along the lateral surface. 

The stresses in a beam of elliptic cross section twisted by constant tangential traction 
are obtained as an application of the theory. 



1. Introduction 

The classical Saint-Venant theory of torsion 
presents a method for the determination of the 
elastic behavior of an orthotropic elastic beam of 
uniform arbitrary cross section twisted by forces 
applied on the end sections, with the lateral surface 
free of stress. In the present paper the theory is 
extended to include a more generally anisotropic 
cylinder (possessing only a single plane of elastic 
symmetry, rather than three mutually orthogonal 
planes of elastic symmetry, as in orthotropic ma- 
terials) built in at one end and twisted by forces 
uniformly distributed along the lateral surface. 

The extension of the Saint-Venant theory of tor- 
sion of cylinders by forces applied on the end sections 
to anisotropic cylinders was made by Lechnitzky 
[3, 4] 3 in 1939; 

The torsion of isotropic cylinders by^forces applied 
on the lateral surface has been considered by Filon 
[1] in 1902, Timpe [11] in 1912, and Zwolinsky and 
Riz [12] in 1939. 

Filon and Timpe showed that the stresses and dis- 
placements in an isotropic circular cylinder twisted 
by forces applied on the lateral surface are derivable 
from a biharmonic stress function. Zwolinsky and 
Riz treated the problem of the isotropic cylinder of 
arbitrary cross section twisted by constant tangential 
forces, and showed that, in this case, the stresses and 
displacements are derivable from two stress func- 
tions, one biharmonic, the other harmonic. For a 
circular cylinder the harmonic stress function van- 
ishes, and the Filon-Timpe solution is obtained as a 
special case of the Zwolinsky-Riz solution. 

The results obtained here are more general than 
those given by Zwolinsky and Riz in two important 
respects. First, the applied forces assumed in the 
present treatment are not necessarily constant or 
tangential to the lateral surface; and second, the 
results are applicable to a wider range of elastic 
materials, including wood, which is orthotropic, and 
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certain crystalline media [5, p. 155 to 160], which 
possess a single plane of elastic symmetry. 

It will be noted that no mention has been made of 
the distribution of stresses over the built-in end of 
the cylinder. This distribution will, in general, 
depend on the mode of clamping. If these stresses 
are prescribed the exact solution of the problem 
presents considerable difficulty because of the com- 
plicated boundary conditions. 

In 1855 Saint-Venant stated as a principle of 
elasticity that two statically equivalent systems of 
forces applied to a portion of the surface of a medium 
will produce essentially the same stresses at point 
of the body not too near the region of application of 
the forces. 

If this principle is accepted, and if it is borne in 
mind that in most applications only the resultant of 
the stresses on the end of the cylinder is known, it is 
seen that a solution of practical usefulness may be 
obtained by ' ^relaxing' ^ the boundary conditions, and 
specifying only the resultant stress on the end section. 

The extra freedom allowed by the relaxation of the 
boundary conditions permits us to obtain a solution 
by the so-called semi-inverse method of Saint- 
Venant, which consists in making certain assump- 
tions about the stress distributions, while still 
leaving enough arbitrariness in the assumed stresses 
to satisfy all the conditions of the problem. 

By the use of this method it is shown that the 
stresses in the twisted cylinder are derivable from 
two stress functions, <t>(Xjy) and yp{x,y), which satisfy 
the differential equations. 
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in the region B corresponding to the interior of the 
cyhnder, and the conditions 



/ 7 T^\ 



M2) 



on the bounding curve C of each cross section. The 
parameters Z, m, 7^, 6^, and hij appearing in these 
equations are functions of the 13 elastic constants of 
the material, N and T are the normal and tangential 
components of the applied surface tractions per unit 
area, v designates the unit normal vector to C, and 
a is given hy 

f {xY'-yX)ds 



In eq (3), X and Y are the x and y components of 
the applied surface traction, and the contour integral 
is the applied twisting moment per unit length. The 
surface integral is proportional to the torsional rigid- 
ity of the cylinder. 

In the case of an isotropic beam the differential 
equations (1), simplify to 



the boundary conditions (2) become 

\dx~V ^^^ ^^' ""^"^(b^^ V ^^^ ^^' """^^^ 

{-^+-)cos(y,v) 



(4) 



(^|__^cos(?/,.) 
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and a is given by 



^=2 



r {xY-yX)ds 



L{''+y'+''4-'-iy' 



(6) 



It is shown, from energy considerations, that the 
boundary-value problems posed above possess a 
unique solution, and several methods of solution are 
discussed. 

As an illustration of the method, the stresses in a 
cylinder of elliptic cross section twisted by constant 
tangential traction are obtained. 



2. General theory 

The first fundamental boundary-value problem of 
elasticity may be stated as follows [8, p. 72]: 

Determine the distribution of stress and the dis- 
placements in the interior of an elastic body in 
equilibrium when the body forces are prescribed and 
the distribution of the forces acting on the surface of 
the body is known. 

The second fundamental problem is similar, except 
that in this case the surface displacements rather than 
the forces acting on the surface are prescribed, and in 
the third, or mixed, boundary-value problem the 
forces are prescribed over a portion of the surface and 
the displacements over the remainder. 

The problem of the cylinder that we are considering 
is of the first kind, where body forces are absent. 

In the following paragraphs we shall review certain 
basic elements of the mathematical theory of elastic- 
ity and present the mathematical formulation of the 
first fundamental boundary-value problem. 

The state of stress in an elastic medium is char- 
acterized by the six components r^x, ryy, r^g, r^g, Tz^j 
Txy, of the stress tensor. The components of stress, 
Xp, Fy, Zr, on an element of surface whose normal 
has the direction V are given in terms of the stress 
tensor by the relations [8, p. 40] 



Yy=T^y cos(x,v) + Tyy cos(y,p) + Ty^ cos(z,v) y (7) 

Z, = rz:c COS(x,v) + Tyz COs(y,v) + Tzz C0s(z,v) . ^ 

The relations (7) may be regarded as defining the 
six components, Txx, Tyy, . . ., T^y, of the stress tensor. 
If we let the vector 7 take on the directions of the 
three coordinate axes in eq (7), we obtain the usual 
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(8); 


X, = r,x, 


■^ Z = TyZf 


Zz^Tzz.^ 















From considerations of static cfiuilihriiim it can be 
shown that 



Y — T^ 7 — Y 



Y^=Z, 



(9) 



and that the defining relations (8) are self-consistent. 
For any material in static equilibrium under a 
prescribed system of surface loading Xy^ Yv, Zp, and 
free of body forces, the stress tensor must satisfy the 
equilibrium equations 



^Txx I ^'^xy I ^'^zr r\ 
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in the interior of the medium, and the conditions (7) 
on the surface, whe^re we now regard X^, Yv, Zp as the 
prescribed surface tractions. 

The strain tensor, e^cxy ^w, ^zz, ^yz, Czx, Oxy, is defined 
in terms of u, v, and w, the x, y, and z components 
of the displacement vector, by the relations 
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It may be shown [5, p. 48] that a necessary and 
sufficient condition that a given tensor may be con- 
sidered as a strain tensor arising from actual physical 
displacements in an elastic material is that the 
tensor satish' the conditions 
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Equations (12) are known as the compatibilit}^ 
equations and are the integrability conditions for 
the system of (11), when the displacements u, v, w 
are to be determined from a given strain tensor. 

The stress and strain tensors in an elastic medium 
are related by the generalized Hooke's law. If the 
elastic constants of the medium are denoted by 
(^a (hj=^,'^j ' ' ' fi) Hooke's law for a medium in 
which the xy-plsme is a plane of elastic symmetry * 
may be written in the form [8, p. 62]: 



^xx = CiiTxx-\- Ci2Tyy-{- CisTzz-\- Ci^Txy 
^yy^ C21'7'tx "T C'22Tyy I <^23'^22 H~ C2Q^xy 
^■zz = CziTxx + Cs2ryy + CisT22 + CsQTxy 



CuTyz-\-C^^Tz_ 



>} 



(13) 



^zx ^ 54'^ yz I (^55'^zx 

^Xy^=C6irxx-{-CQ2Tyy-\-CQsTzz-\-CQQTxy J 



where Ci2==C2i, ^13=^31, Cie^Cgi, 023=032, 026=^62, 

0.36 "=063, 045 = 054. 

For isotropic materials, the number of independ- 
ent elastic constants reduces to two, and Hookers 
law has the simpler form 



^xx TpTxx Tfi\'^yy\'^zz)j ^ijy — rp'^yy' 



TpKjzzKTxxjf ^zz — Tp '^ ^ 
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where E is Young's modulus and o- is Poisson's ratio 
for the material. 

We may now reformulate the first fundamental 
boundar3^-value problem of elasticity in mathemat- 
ical terms: 

Determine a set of six functions, the components 
of the stress tensor, which satisfies the equations (10) 
in the interior of the elastic medium and the relations 
(7) on the surface of the medium, where Xy, F^, and 
Zv are prescribed functions, and such that the strain 
tensor calculated from these functions by the use of 
(13) satisfies the compatibility conditions (12). 



3. Formulation of the problem 

Let the coordinate system be chosen as shown in 
figure 1, with the built-in end of the cylinder in the 



< The most general elastic medium is characterized by 21 independent elastic 
constants. Here we consider a medium whose elastic properties are identical in 
any two directions symmetric with respect to the ry-plane. This type of elastic 
symmetry reduces the number of independent elastic constants to 13. 
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jplane z=0, and the free end in the plane z=h. The 
boundary curve of any cross section will be denoted 
by C; the region enclosed by C and the unit normal 

to C will be denoted by R and v. 
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The prescribed surface traction T (force per unit 
area) is assumed to be independent of z, directed 
parallel to the (x,y) plane, and such that the result- 
ant on any cross section is a twisting moment of 
magnitude M (torque per unit length). 

This condition requires that the components X 
and Y of the applied surface traction satisfy the 
relations 

r Xds=0, r Yds=Q, r (xY-yX)ds=M. 

(14) 

In order that the first fundamental boundary 
value problem of elasticity may be formulated in 
such a way that the solution is unique, it is necessary 
that the stresses be prescribed over the entire sur- 
face, subject, of course, to the conditions of static 
equilibrium. This would require that Zx and Z^, the 
components of the clamping stresses on the end 
2=0, be prescribed, subject to the conditions 

f ZJS=0, f ZydS=0, ( (xZ,-yZ,)dS = Mh. 

JR JR Jr 

(15) 



In this case the boundar^^ conditions on the end 
sections would be of the form: 



(a) on the end s=0: 

(b) on the end z=h: 



^ zz — ^x^ 



.=0, 



Tyz — ^y\ 



r = 0. 



(16) 



Inasmuch as Ti^ and Z^ are unknown (and depend 
upon the mode of clamping), we shall replace the 
exact boundary conditions (16) by the relaxed 
conditions. 



(a) on the end 2=0: 

TzzdS= XT2zdS= yTzzdS = 
Jr Jr Jr 

I TzJS= \ TyJS = 

Jr Jr 

{xTzy—yT2x)dS = Mh; 
Jr 

(b) on the end z=h: 

TzzdS= xr2zdS= yTzzdS = 
Jr Jr Jr 



r = 



z = 0. 



>■ (16a) 



These are conditions on the resultant stresses, rather ^ 
than on the stress distributions. 

The boundary conditions on the lateral surface, 
that is, on the curve C, are 



Txx cos {x,v) + Txy COS (yyv)=X' 
Txy COS {x,v)+ryy COS {yjv) = Y 
Tzx COS (x,v) + Tyz cos (y^v) = 0.. 



(17) 



Our problem, then, is to determine a system of 
stresses Txx, Tyyy . . , j Txy which satisfies the equilibrium 
eq (10), the relaxed boundary conditions (16a) and 
(17), and for which the strain system calculated by 
the use of Hooke's law (13) satisfies the compatibility 
equation (12). 

The solution of the problem as thus formulated is, 
of course, not unique, but, according to the principle 
of Saint-Venant [8, p. 95], the stresses given by any 
two solutions will be essentially the same at points 
of the cylinder not too near the end sections. 

We shall make use of the Saint-Venant semi-inverse 
method of solution [8, p. 100]. This consists in 
making certain assumptions about the stress distri- 
bution, and still leaving sufficient freedom to permit 
the satisfaction of the differential equations and the 
boundary conditions. If these assumptions are 
carefully made, a unique solution (unique under 
these assumptions) exists. 

It is apparent that the torsional couple trans- 
mitted across each cross section is a linear function 
of the distance from that cross section to the free 
end of the cylinder. Hence, it is natural to assume 
that the stresses in the cylinder are given, in part at 
least, by a Saint-Venant torsion function. 

There is a type of elastic deformation known as 
plane strain [7, p. 231], in which the displacements , 
are of the form 

u=u(xyy), v=v(xjy)^ w=0. 
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The stresses in this type of deformation are derivable 
from an Airy stress function and can be chosen to 
satisfy very general boundary conditions of the type 
given by (17). 

We shall assume that the stress system in the 
twisted cylinder may be derived, in part, from a 
superposition of the torsion, solution and the plane 
strain type solution. 

For a deformation to be of plane strain type, it is 
necessary that Tzz^^O. Because a variable tension 
on the end sections may produce a bending or exten- 
sion of the cylinder, we shall further superpose the 
bending and tension type solutions. 

Thus, if we denote by rtj any one of the six com- 
ponents of the stress tensor, we m.ay write 



TiJ 



= ^(1)J_^(2) I (3) 



(18) 



where a is a constant to be determined by conditions 
of static equilibrium, r/J^ is derived from an Airy 
stress function, r/f is derived from the Saint-Venant 
torsion function, and r^f is required to satisfy the 
condition that the end surfaces be free of stresses 
that would produce bending moments and extensions. 
The stress system r^f may be written [7, p. 231] 



dxdy 



.(!)__ 



5^ 

57/ 



. (1) ^ Q 

xz ^ ) 



r(l) = < 



(19) 






where i^[ = -^{x,y)] is the Airy stress function and ^, 
q, and r, are functions of the elastic constants. We 
note that the stress system r/f satisfies the equilib- 
rium equation (10) identically. 

The stress system r/f has the form [8, p. 101, 108] 

(20) 
where -4, B, (7 are determined by the three relations 

J^ r^^dS=0, ~J^ XT^zdS=Q, j^yT,,dS=0, 

(21) 

which state that the resultant bending couples and 
tension vanish on the end sections. 

The system r/f is a modification of the Saint- 
Venant stress system r-,-. Sokolnikoff [8, p. 214] 
gives the stresses r^'y in an orthotropic beam twisted 
by end couples as 



'■■-=" (sJ+')+'"(s|-'') 



T ^^^ T ^^^ T 



ZZ—Tl,= Oy 



(22) 



where I, m, n are functions of the elastic constants 
and (t>[ = (t>{Xjy)] is the torsion function. Because 
we wish the torsional couple to vanish on the face 
z=h, we shall set 






(23) 



The third of the equilibrium equations (10) will be 
satisfied if <j)(x,y) satisfies the condition 



dV 
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dx^ bxby ' '"5?/^ 






(24) 



In order that the first two equilibrium equations 
in the set (10) be satisfied identically, we set 



T«'=^.^+fa;^ 



--n4>- — y\ 



■if=^x'-l:y'+ni<i>. (25) 



To simplify the satisfaction of the compatibility 
conditions (12), we choose 



rf: = -<l>+^^-ly\ 



(26) 



with s and t as yet undetermined. 

Our assumed stress system, then, is the following: 

1 7 . c)V_i_^ 2 

1 , dV ^ 2 
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r (27) 



267 



with the constants a, /, m, n, ^, g, r, s, and i as yet 
undetermined. 

We next determine ]), q, r, s, and ^, in such a way 
that Czz will be independent of i^(x,?/) and of the 
terms in x^ and y^. The insertion of the stresses 
given by (27) into the third of (13) leads to the 
following expression for 6,2. 



1 



c,^{A+Bx+Cy) + 
1 



(C3im+C33s + C36ri)x2— 

(C32m + C33^ + c^J)y^ + (C332? — C32) 






2+ 



(Csgg + Cae) 



dV 



"(^33^ ^^^^d^/^' 



dxdi/ 
If we set 

s= — (csim+CsQn)/cs3, t= — (^32^+^36^/^33, 



(28) 



(29) 



the last five terms in the right-hand member of (28) 
will vanish. 

If the values for p, ^, r, s, and t given in (29) are 
substituted in the stress equations (27), and the 
strains are then calculated from Hookers law (13), 
we have the following expressions for the components 
of the strain tensor: 
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bV 



dV 



dV 
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C23(^+£:r+Gy) 



C44^+C45m)('^+a; )+ 
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(30) 
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(30) 
Con. 



where 



On — bii — 



CijCiz 
e3j(^ZS 



C33, ii= 1,2,6 



(31) 



6i=CnHci2^— c^3+c^6m, i= 1,2,3,6. 



If the components of the train tensor given by (30) ^ 
are inserted in the compatibility equations (12), it 
is seen that all but the fourth of the compatibility 
equations are satisfied by the assumed stresses if l^ 
m, and n can be chosen to satisfy the six equations 



Cun— €5^1=0 

^45^ + ^55^ = 

^(C55+C3l)+m(2C45 + C36)+n(C44 + C32)=C33 
/(C3l)+m(2C45 + C36)+n(2C44 + C32)=C33 
^(2C55 + C3l) +m{2C45 + Cz6)+n(C32) =Czz. ^ 



(32) 



(It is easily seen that only three of these six equa- 
tions are linearly independent and that the compati- 
bility equations are satisfied if we set 

l=Czzcj0j m=—C3scJ^y n=czzc-j^, (33) 
where 

/3 = C44(C55 + C3l) +^55(044 + ^32) — C45(2C45 + C36). 

For convenience in calculation we rewrite the ex- 
pressions for s and t given by (29) in terms of the 
elastic constants, using the relations (33) : 

S= (C31C45 — C36C55)//?, t= {Cz2Cio — CzQCu)/^. (34) 



The fourth of the compatibility equations (12), 
the only one that is dependent on \l/{x,y), requires 
that \p{x,y) satisfy the relation 
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dV 
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by' 



b'<l> 
bxby 



bi^,+ h,n-bj. (35) 



In the next section we shall consider the boundary 
conditions on the functions (l)(x,y) and ypix^y). 
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4. Boundary conditions 

It has been shown in the previous section that the 
system of stresses given by 






Txx=i<t> — ^+'-K- X^ 
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will satisfy the equilibrium equations (10) if <i>{x,y) 
satisfies 



dV 



^H , c)'<t> 






(37) 



and will satisfy the compatibility conditions (12), 
provided that 

13 = Cu (C55 + C31) + C55 (^44 + C32) — C45 (2^45 + Csq) , 

and \ly{Xjy) satisfies the differential equation 



(38) 
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ICz: 



^,i=l,2,6, 



and 



hi=Ciil+Ci2n—Ci^+CiQm, 'i== 1, 2, 3, 6. 



(40) 



Since both of the functions <t>{Xjy) and yp{x,y) 
appear in (39), it is apparent that we must first 
determine ^{x,y). Fortunately the third boundary 
condition of the set (17) is independent of yp(x,y), 
so that the two boundary value problems are 
separable. 

If the values of Tyz and Tzx given by (36) are 
inserted into the third of (17), we see that this 
boundary condition will be satisfied if 

(ly —mx) cos {x , v)-\-(my—nx) cos (y, v) (41) 

on C It is shown by Sokolnikoff [8, p. 215] that the 
change of variables 



m 



^=-yjln—m'^x, v = y — fX 



^,(^^^)_>c:m!^(j^^) (42) 



I 

will reduce the systems (37) and (41) to the simpler 
system 

y</>Ur/)L^V(«,^) 



d^ 



dn' 



= 



(43) 



d ^r 



di 



0'(^> 11) = ^ cos (^jv)—^ COS (tj,v)=0 on C\ 



The problem of determining <l>{x,y) is thus reducible 
to a problem of Neumann. Since the right hand 
member of the second of equations (43) is a perfect 
differential, the determination of (t>^ {^,r() is possible, 
and is unique to within an additive constant. 

It should be pointed out that the quantity 
(In—m^) is positive for all elastic media, since from 
the defining relations (33), we observe that 



In m ="32" (^44^56 ^45) 



(44) 



and the positiveness of (C44C45— C45) is guaranteed by 
the positive definite nature of the strain energy form.^ 
Once the function <l){x,y) has been determined the 
right-hand member of (39) is a known function. 

Let us express the surface tractions X and Y in 
terms of the tangential traction T and the normal 
traction N, by the relations 

X=N cos {x, v) — T cos (t/, v) , 

Y=N cos {y, v) + T cos {x, v) , (45) 



5 See appendix. 
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If we make use of the relations (45), the stresses 
given by (36), and the relations 



cos (X, v)=j^y 



cos{y,v)=-j^> (46) 



M47) 



the first two equations of (17) may be rewritten in 
the form 

(^x''—^y^+m^+-Jcos{y, v), 

(rK^-^y^-— )cos(2/, v), 

on C. Let us denote the right-hand members of 
these equations by F2ix,y) and Fi{x,y), respectively. 
If we can find a particular (single-valued) solution, 
^*{x,y) of (39), we must then solve the homogeneous 
system: 



da;*" 



dV 



5V 



^22 ^Ti- 3 ^26 s:3k:::+2(Si2+ hj ^^25^2' 



dx^dy 



36, 



bxby 



+ 611 






ds 



Y (48) 






It is easily deduced from the results given by 
Sokolnikof [7, sections 66 and 95] in his discussion of 
the theory of anisotropic plates that the system (48) 
possesses an essentially unique, single-valued solution 
f (a:, y)j provided that (a) the roots of the character- 
istic equation 

622X^-3626X^+2(612+666) X'-36i6X+6ii-0 (49) 

are all complex, and (b) 

r F2ds = 0, r Fids = 0, r (xFi-yF2)ds = 0. 

(50) 

It was shown by Lechnitzky^ [3] and [4], through 
energy considerations that the roots of eq (49) are 
always complex. Thus condition (a) is fulfilled. 

We shall test the functions Fi(x, y) and F2{Xf y), 
in the first two equations of (50). These equations 
become 



« Lechnitzky considered a similar equation, arising in the study of anisotropic 
plates, whose roots are the negatives of the roots of eq (49). 



=- [N cos (x, v)—T cos {y, v) ]ds, 
oijc 

y (51) 

Jc(S ^'""2 ^'+^ V ^^~(^*~Y ^ / ^^ 

=- N co^{y^v)-\~T Q,o^{x,v)ds. 
oiJc J 

The right-hand members, by the use of eq (45) are 
seen to represent the components of the resultant 
force per unit length. 

By making use of Green's theorem, 

we rewrite the left-hand members of (51) in the forms 



LH%-'><»V'- 



(53) 



Since 4>{x,y) satisfies the differential equation (37) 
we may rewrite the expressions (53) in the forms 



(64) 



da; 



Once more we use Green's theorem (52), this time on 
the expressions (54). These now have the forms 



ds 



J ^1 1^ 5^+*^ 5^-^2/ + ™^| cos(a;,j')+ 

|m ^+71 ^-mi/+wx| cos {y, v)^ 

L^W^ b! +^ |f-^2/+mx| cos {X, .)+ 

\m ^+^ ^— my+nx) cos {x,v) ds 



(55) 
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The expressions in brackets vanish, by (41), and 
hence the first two of our integrabihty equation 
(50) become, in this case: 

r Xds=0, r Yds = 0. 

These are the conditions which the prescribed 
components X and Y of the appHed stress were 
assumed to satisfy in section 3, hence the first two 
integrabihty conditions (50) are satisfied. 

Let us now consider the third equation of (50). 
We express N and T in terms of X and Y and see 
that we must have 

-( {xY-yX)ds 

We apply Greenes theorem to the right-hand member 
of (56), to obtain 

- f (xY-yX)ds 

Now, the left-hand member is proportional to the 
twisting moment per unit length, M, that is, 

M= r {xY-yX)ds, (58) 

Let us denote by D^ the expression ^ 

Z?'= nx^ — 2mxy-{-ly^+x ( m -zr — \-n r— j— 

Then, if a is taken as 

a=MID\ (60) 

the third equation of (50) is satisfied, and a single- 
valued solution of (48) exists. 

We have satisfied the boundary conditions on the 
lateral surface; next we shall consider the relaxed 
boundary conditions on the end sections. 

On both end sections we require that there be no 
resultant forces. Hence we require that 



^ The quantity D' is proportional to the Saint- Venant coefficient of torsional 
rigidity i>, which arises in the theory of torsion by end couples. 



f r.,rfS = 0, f T^^dS=0, f T.,rfS = 0. 

JR JR JR 

The first of these conditions can be satisfied by 
adjusting the parameters A, B, and C. The second 
and third are easily shown to hold by considering 
(53) through (55). 
On both ends we require that the bending couples, 

z z 

Mx and M^, vanish. These requirements can be 
met by adjusting the parameters A, B and C, so that 



-X 



M,= yT,,dS=0, 



JR 



xt^JS = 0. (61) 



On the end 2=0, we require that M^^M-h and on 

z 

the end z=h, we must haveTM2=0. The second 
condition is easily seen to be satisfied since rxz= 
Tyz=0 for 2=\i. 

z 

Since Mz for ^=0 is given by 

M...„ £{.[»(!+.)+» (II-,,)]- 



dS, (62) 



and the integrand is equal to the integrand of (57), 
it follows that 



Mz = haD' = hM. 



(63) 



Thus the stresses given by (36) can be made to 
satisfy the boundary conditions (16a) and (17), the 
conditions of static equilibrium (10) and the com- 
patibility conditions (12) by proper choice of the 
constants I, m, and n, and the functions (l)(x,y) 
andiA(x,?/). 

In the next section we shall specialize the equations 
to the isotropic case, and in section 6 we shall 
illustrate the use of the method by calculating the 
stresses in an elliptic cylinder twisted by constant 
tangential traction. 

5. Isotropic beam 

If the beam is isotropic the differential equations 
and boundary conditions are considerably simplified. 
In this case the elastic constants Cij may be expressed 
in terms of E, the Young's modulus, and a, the 
Poisson ratio, as follows. 



Cii — C22 — C33 — 



E 



C23 — C32 — ^21 — C12 — Ci3 — C31 — — -^ 



C44 — C55 — Cee — 



l+<r 



y (64) 



E 



^16 ^61 <^26 Cro C 



-Cqz C45 — C5 
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The elastic constants I, m, n, hi, bfj are given by 
m=0 

h^={l + a)IE 

b,,= -ia+a')/E 

^16 = ^26=0 

b,=b,=-b,= {l + o)/E 
66=0. 
The function <t>ix,y) must satisfy tlie conditions 

V'<t>(x,y)=0 



(65) 



-p=2/ cos {x,v)—x cos (y, v), on C\ 
and lACa!,!/) must satisfy 

v'Hx,y)=o 

(67) 

The torsional rigidity D^ is calculated from the 
defining equation 

and is related to the Saint-Venant coefficient of 
torsional rigidity D by the relation 

D'=D(l + ^)IE. (69) 

The stress system in this case has the simpler form, 



1 _±_^ 






2 cix^ 



1 



i r,,= -cV^-<i>+{A+Bx+Cy) 



^r,.=i(A-.)(g+a;) 



(70) 



1 1 . 2 2N , t>V 



(70) 
Con. 



1 

-e 
a 



and the strains are given by 

{A+Bx+Cy) 

z^ii, 111 ujo Uj uy Jjj 

1 



'2£''^"^ £■ dx^ E by' 



2E 



<t>- 



■—^<l>+^iA+Bx+Cy) 



^e..=^(/.-^)(^--2/^ 



Y^-, 



.i±^^V+l±^(^2_^2) 



6. Stresses in an elliptic cylinder 



(71) 



As an illustration of the theory developed in the 
preceding sections, we shall calculate the stresses in 
an anisotropic cylinder of elliptic cross section 
twisted by a constant tangential traction T. 

Let the equation of the elliptic boundary, (7, of the 
cross section of the cylinder be given in the paramet- 
ric form 

x = acos 6, y=b sin 6, (72) 

Then, 

cos {x, v) : cos (i/, v)= b cos 6: a sin 6. (73) 



Let us set 



and 



<t>i^,y)= S Aj^'y'y 

i-hj<2 



ii=0,l,2. 



(74) 



i(x,y)= S B.jx'y^, i,j= 0,1,2,3,4, (75) 



and seek to determine the coefficients Atj and Bfj in 

such a way that (37), (39), (41) and (47) are satisfied. 

We first calculate Mfrom the defining relation {5S), 



M= i {xY-yX)ds. 



(76) 



If we insert in (76) the expressions for FandX given 
by (45), we find that M is given by ^ 

8 Although N is assumed to be zero in this discussion, we include it in eq (77) 
to show that a constant normal traction produces no twisting moment. 
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M=^J^{'-yN+xT) cos (j, v) + {yT+xN) cos [y, v)]ds, 

(77) 

By making use of Green's theorem (52), we may re- 
write (77) in the form 



Jr 



dS=2'irabT, 



(78) 



We insert the expression (74) into (37) and (41) 
and, after some manipulation obtain 

(79) 

as the solution of the systems (37) and (41). 

Using the value of (f){x, y) given by (79) in the defin- 
ing equation (59) we find that D' is given by 



j^f _'Ka'^h\nl—m'^) 
~ ihH+a'n) 

Hence, from (60), it follows that 

J2.{hH+aH)T 

We next insert the expression (75) for i/(x,y) into 
the two equations (47) and after setting x=a cos d, 



(80) 



(81) 



y=b sin d, and making use of (73), \ve\)btain a pair 
of polynomials in sin 6 and cos 6. Tliese are cwpanded 
in Fourier series and the coefficients of the various 
terms are equated to zero. We thus find that, in 
order to satisfy the boundary conditions, \f^(x,y) 
must be given by the expression 

24a'b'(a'n+bH)'x(^(x,y)=^b'CSlmb'-5mna'+y)x^+ 

4a'b'(b'm'-a'n')7?y+2a%^yxY 

+4a'b'(b'P'-a'm^)xy'+a\5lmb'-Smna'+y)y'- 
2a'b^(3lmb'-37nna'+y)x'- 

2a'b\Slmb'-3mna^+y)y\ (82) 

where 7 is a constant, as yet arbitrary. 

If the value of \p{x,y) given by (82) is inserted in 
the differential equation (39), it is found that 7 must 
satisfy the relation 

{3b22b'+2a'b'(bi2+b,,)+3br,a'}^y 
=3b22b\bmna'-3lmb'')+^b2^d'b\¥7Yi'-aV) + 
%u(i''b\bH''-a'm'') + 3bna\3mna'-blm¥)- 
6b2a'b'm+6b,a'b'(a'n-bH)+6bia'b'm+ 

3a'b\a'n+bH) (b2,n-bj). (83) 

Inserting the values of 0(x,7/), a, and yp(x,y) given 
above into (36) we obtain the stress system 



r,,^ (37nna'-3lmb'-y)(bV+3a'y^'-a'b^)-Qa^\ln-m')xy 
T 3a'b\ln-'7n^) 

ryy_ {3mna'-3lmb'-y){3bV+aY-a'b') + Qa^\ln-m'')xy 
T" 3a'b\ln-m^) 

r^^ 3(Zn-m^)(6V-aV)-2(3mna^-3Zm 6^-7)0;?/ 
T~ 3a^b\ln-m^) 

Tyt_ 4:X{h—z) 

rzx_ 4:y(h — 2) 
T~ b' 

^=v'x'+q'xy+r'y'+A+Bx + Cy, 
where 
^_3b^Cz2(3lmb^—5mna^+y)—3a^CsQ(b^m^—aV) + a^Cziy Csim+CzeU. b^m 



V 



Qa%,{bH+a'n) 



^33 



(bH+a'n) 



I ,^ 3c^2{b'm'-aV)- c,,y+3c^i(bH''-a'm') + 3c,,(a'n-bH) 



3c^,(bH+a^n) 
,_b^Cz2y-—3b^c^Q{bW—a^m^) + 3a^Csi{5lmb^--3mna^+y)^Cz2''n+CzJ^ 



Qb%,(bH+a'n) 



C33 bH-\-a^n 



(84) 



(85) 



273 



To evaluate the parameters A, B, C, in order 
that the conditions 

T^^dS={), I Xr^^dS=Q, yT,,dS=0 

Jr jr Jr 

are satisfied we write 

f {^'x^+(i'xy+r'y^-{-A+Bx+Cy)dxdy=Q 
Jr 

I (p^x^-]-q^x^y-^yxy^-{-Ax+Bx^-}-Cxy)dxdy = 

j {p^x^y+q^xy^+r^y^-\-Ay+Bxy+Cy^)dxdy=0. 
JR -J 

(86) 

If we denote by R, Ix, and ly the area and the 
moments of inertia of the elhptic section, the system 
(86) can be written in the form 

p'Iy+r'I,+AR=0, Bly=0, CI,=0, (87) 

For the ellipse, 

E=7rao, I^=—j^7 Iy=—j — (88) 

Hence the parameters A, B, and C are given by 

^__ay + 6V^ 5=0, C=0, (89) 

and the proposed problem has been completely solved. 

For the isotropic elliptic cylinder the solution may 
be simplified by using the relations given in section 5. 

The solution for this case is given by: 

<t>{x,y)= ^2j^^2 ^y 



D'= 



ira'P 






a'¥ 



7 = 



Txx_ 2xy 



2xy 



Txp_b'x^—a^y^ 



(90) 



Ty^ 4x(h — 2) 




T a^ 

Tzx 4y{h—z) 


(90) 


T b^ 

Tzz_2cr(a^—b^)xy 


Con. 



a'b' 



For an isotropic circular cylinder of radius a the 
results are given by 



<l>ix,y)=0 

4 

8T 



'/^(^, y) 



xy^—x^y 



yj- 


24 


Txx 


2xy 


T 


a? 


'^m_ 


2xy 


T 


a' 


Txv 


x'-y' 


T 


a? 


ryz 


Axih—z) 


T 


a" 


Tzx_ 


4:y(h-z) 



(91) 



— =0 



7. Methods for obtaining solutions of the 
boundary-value problems 

The determination of the stresses and displace- 
ments in the twisted cylinder has been made to 
depend upon the solution of the two linear^partial 
differential equations. 






dx' 



bxdy dy^ 



(92) 



and 



t>^^ -5^4-3626 5^35^+2(612+6..) 5^25^-36.. 5^5p+ 
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together with the appropriate boundary conditions, 
(41) and (47). 

It was shown in section 4 that (92) can, by a simple 
change of variables, be transformed into the form 



V'<I>{X,7J)=0, 



(94) 



and the various methods of potential theory are avail- 
able for the solution of this problem of Neumann. 
An elegant method of solution of the Neumann 
problem, due to MuscheliSvili, is discussed in 
Sokolnikoff [8, p. 170]. 

The determination of the function ^{x,y) is con- 
siderably more difficult. An equation similar to (93) 
arises in the theory of anisotropic plates and a great 
deal of work has been done on the problem of solving 
this equation. A series of papers by Green and 
Taylor [2] gives a method of solution of the system 
(93) -(47) for the special case (corresponding to the 
ortho tropic beam) , where 



b26 = bl6 = 0, 



(95) 



A perturbation method of solving the system (93)- 
(47), which consists in expressing the solution as a 
series of biharmonic functions, is described by 
Sokolnikoff [9]. Another, more elegant method, 
essentially due to Muscheli^vili, depending on con- 
formal mapping is also described bv Sokolnikoff 
[7, p. 375] and [10]. 

If the beam is isotropic, (93) has the simple form 

VV(^,2/)=0, (96) 

and the solution is easily obtained in the form 



Hx,y)=R[m + 2g(2)l 



(97) 



where /( 2) and g{z) are analytic functions of the com- 
plex variable 2:[ = x+%], which are chosen so that the 
function i/{x,y) given by eq (97) will satisfy the 
boundary conditions (47) [10]. 

8. Summary 

The determination of the state of stress and the 
displacements in a built-in anisotropic cylinder 
twisted by forces uniformly applied along the lateral 
surface has been made to depend on the solution of 
two boundary-value problems. 

The existence of the solutions of the boundary- 
value problems has been demonstrated, and several 
methods of solution of the boundary-value problems 
are indicated. 

As an illustration of the method, the stresses in an 
anistropic beam of elliptic cross section twisted by 
constant tangential traction were obtained by assum- 
ing that the stress functions were second- and fourth- 
degree polynomials. 



In a succeeding paper it is planned to extend the 
theory to the case where the applied surface tractions 
vary along the length of the cylinder, that is, where 
the boundary conditions are given in the form.^ 



Y=g(z)Y{s) on a 

z=o. 



(98) 



It is hoped that the theory can be further extended 
to the problems of extension, bending, and flexure of 
anisotropic cylindrical beams by forces distributed 
on the lateral surf ace. ^^ 

The boundary conditions for these problems are: 

(a) Extension: 

X=0 on C 

r=0 on C, 

and Z is prescribed on C subject to the conditions 

f Zds9^0, \ xZds = 0, \ yZds = 0, (99) 

%J C u C *) C 

(b) Bending: 

Z=0 on C 

F=0 onC, 

and Z is prescribed on C subject to the conditions 

r Zc/s=0, r xZds9^0, r ^/Zc^s^^O. (100) 

(c) Flexure: 

Z=0 on C, 

and X and Y are prescribed on C subject to the 
conditions 

{ lXds9^{), r Ydsy^{), r {xY-yX)ds=0, 

(101) 

It is apparent that a superposition of solutions 
corresponding to the boundary conditions (98) 
through (101) can be made to satisfy quite general 
boundary conditions on the lateral surface. 



9 Zwolinsky and Riz (12] treated the case (for isotropic cylinders) where f{z)—z, 
N=0, T= constant on C. 

10 The problems of extension, bending, and flexure of isotropic (and orthotropic) 
cylinders by forces applied on the end sections have been solved by Saint- Venant 
[8 p. 100]. 
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10. Appendix. Strain Energy and 
Elastic Symmetry 

The strain energv, W, of a strained elastic medium is 
defined by [8, p. 88]^ 

2 W= T:cXeXX -h TyyByy ^ T.^S.^ 4" ^Ty^By;, 4" ^T.^B,^ + 2Txyey^. (102) 

In the linear theory of elasticity the strain energy is 
assumed to be given also by a positive definite quadratic 
form in the components of the strain tensor by 

W = i^n^L ^-hnexxeyy-\- feis^xz^zz + buexxeyz + huexxCgx + bi^Cxxexu 

+ 2 &22e^y 4- h2Ze yyBgz + h2\eyyeyg + 625^yi,62X + &26^I/y^Xl/ 

+ h^zzel + 634^2,6^2 + hzf^ezzCgx + hueg^exy 

+ ibuelt+bi^eygezx -\-bi^eyzexy 

-{-hbbhe^-^-h^egxexy 

(103) 



The positive definite nature of the form (103) imposes certain 
requirements on the signs of the elements of the matrix 



pu 


612 


613 


6.4 


615 


&.6^ 


612 


&22 


&23 


624 


625 


626 


613 


623 


633 


634 


635 


636 


6h 


624 


634 


644 


645 


646 


fc.5 


&25 


635 


645 


655 


&56 


J>\t, 


h. 


&36 


646 


6» 


hm-^ 



(104) 



An elastic medium is said to possess elastic symmetr\^ in a 
plane if a reflection of coordinates in that plane leaves the 
strain energy form (103) invariant. Let us consider elastic 
symmetry in the a;?/-plane. The form (103) must be invariant 
under the transformation 



y'=y 

z'= —z. 



(106) 



The components of the strain tensor transform according 
to the relations 






^z'z' — ^zi 



Cz'x' — ^zi 



(106) 



and if the form (103) is to be invariant under the transforma- 
tion, we must have 



&14 = &15 = &24 = &2o = ^35 = ^46 = &56 = 0. 

The matrix (104) then has the form 



(107) 



rhn 


hn 


6l3 








fel6">| 


hu 


h2 


&23 








?>26 


&13 


^23 


^33 








^36 











644 


&45 














645 


655 





<bu 


b2. 


636 








b&6.^ 



(104a) 



Since the matrix (104a) is obtained from a positive definite 
quadratic form it follows that 



644 645 

bio L>55 



>0. 



(108) 



By differentiating (102) and (103) with respect to Cxx^ and 
comparing, Hooke's law is obtained, expressing the stresses 
in terms of the strains. 

If this expression of Hookers law is solved for the strains 
in terms of stresses and the coefficients are denoted by 
cn, . . •, CfiB as in section 2, it will be seen that the determinant 



C44 C45 
C45 C55 



(108a) 



is a positive multiple of the determinant in (108), and hence 
is positive, as asserted in section 4. 



Los Angeles, May 26, 1951. 
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